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ABSTRACT 

Although the two moons of Mars, Phobos and Deimos, have long been thought to be captured 
asteroids, recent observations of their compositions and orbits suggest that they may have formed 
from debris generated by one or more giant impacts of bodies with 0.01 x target mass. Recent 
studies have both analytically estimated debris produced by giant impacts on Mars and numerically 
examined the evolution of circum-Mars debris disks. We perform a numerical study (Smoothed Particle 
Hydrodynamics simulation) of debris retention from giant impacts onto Mars, particularly in relation 
to a Borealis-scale giant impact (E ~ 3 x 10^® J) capable of producing the Borealis basin. We find that 
a Borealis-scale impact is capable of producing a disk of mass ^ 5 x 10^° kg 1 - 4 % of the impactor 
mass), sufficient debris to form at least one of the martian moons according to recent numerical studies 
of martian debris disk evolution. While a Borealis-scale impact may generate sufficient debris to form 
both Phobos and Deimos, further studies of the debris disk evolution are necessary. Our results can 
serve as inputs for future studies of martian debris disk evolution. 

Subject headings: Mars; Mars, satellites; Impact processes; Satellites, formation 


1. INTRODUCTION 

The origin of the two moons of Mars, Phobos and 
Deimos, is widely debated. One prevailing theory 
is that the two moons are captured asteroids, as in¬ 
dicated by their size , shape, and prior observations 
of their coinpositi on (IMurchie et al\ II99H iBurnsI 119781 
Forget et al.ll2008ll. However, recent obse rvations (e.g., 
Giuranna et aLl 120111 : iWitasse et al\ 1201411 suggest that 
Phobos and Deimos may have formed from debris gen¬ 
erated from re-accretion after one or more giant im¬ 
pacts of bodies with ^ 0.01 x target mass, a process 
that could more readily explain the composition, den - 
sity, and orbits of the martian moons (|Rosenblattll2011ll . 
Compositionally, observations from both the European 
Space Agency’s Mars Express mission and NASA’s 
Mars Global Surveyor mission show that Phobos and 
Deimos have thermal infrared spectra in poor agreement 
with the primitive carbon-rich materia ls normally as¬ 
socia t ed with carbonaceo us chondrites (iGiuranna et al.\ 
1201 H IWitasse et ^1201411 . Additionally, minerals were 
detected on Phobos that are also present on the mar - 
tian surface (jGiuranna et al.l[201lUWitasse et alll2014ll . 
suggesting that the material that formed the martian 
satellites could have been ejected from Mars. Further¬ 
more, the Mars Radio Science Experiment on board Mars 
Express mission measured the density of Phobos to be 
1.87 g cm~^ with a porosity of 25-35% I Andert et al\ 
l201Ctl^tzold et a7ll2C) 14 [Rosenblatt et aZ.ll2010[l — such 
low density and high porosity are uncommon for most 
asteroids. An asteroid with such a composition would 
not be expecte d to remain co herent if it was dynami¬ 
cally captured (lGraddocldll994ll . although it could be ex- 
pected to form by mater ial re-accreted from a debris disk 
(|Richardson et aDl2002f) . Re-accretion could also explain 
the circularity of the orbits of Phobos and Deimos; if the 
martian moons were captured objects, they would be ex¬ 
pected to have more eccentric orbits and randomly ori¬ 


ented orbital planes - and although the orbit of Phobos 
could be circularized via tidal perturbations from Mars, 
Deimos is too far away to be strongly affected by such 
perturbations. Additionally, moons that formed from 
the same impact that produced the current martian spin 
would be expected to or bit near the equatorial plane, a s 
Phobos and Deimos do (|Rosenblatt and Charnoz l l2012fl . 

The observations discussed above makes a strong case 
that Phobos and Deimos formed via re-accretion, but 
this process requires the ejection of sufficient satellite¬ 
forming material into circum-Ma rs orbit, most likely via 
a giant impact (jGraddockl 1201111 . Giant impacts were 
a common process in the late stages of planetary for¬ 
mation, and several large basins between 3300 and 4500 
km in diameter suggest impacts that could have ejected 
significant debris into martian orbit. The largest pos¬ 
sible impact site is the 7700 km-diameter Borealis basin 
that encompasses almost the entire northern hemisphere. 
Although whether the Borealis basin formed from exo¬ 
genic or endogenic processes is still debated, recent anal¬ 
ysis of the martian crust suggests that a giant impact 
could explain the elliptical shape of the basin and the 
sharp crustal thickness grad ient on the basin perimeter 
(I Andrews-Hanna et aLl[2008ll . In addition, two indepen¬ 
dent studies numerically simulated Borealis-scale giant 
impacts (E ^ 1 — 3 x 10^® J) and found that such an im- 
pact could have formed a structure simila r to the basin 
(iMarinova et a/.ll2008t iNimmo et a/.ll200^ . 

A Borealis-scale impacto r would likely have ha d a mass 
of ^ 0.026 Mars masses (jMarinova et al\ |2008|1 . allow¬ 
ing the ejection of substantial material into martian or¬ 
bit. The total mass of ejected material necessary to form 
Phobos and Deimos is unclear, and is complicated by the 
large number of elongate craters on the surface of Mars 
that could have formed via impacts by add i tional short¬ 
lived martian satellites (|Schultz et a/.lll98l : iStrom et al\ 
Il992f) . The total mass of past and present martian satel- 
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lites, Msats, is estimated to be between 8.3 x 10^^ and 
1.5 X 10^® kg lx 10“® — 2 X 10“® Mars masses), 
based on the appli cation of the Gault and 7r-group 
crate r scaling laws (jGaultl 11974 iGanlt and Wedekindl 
Il977t IMeloshlll989ll to various estimates of the number 
of elongate martian crate rs attributed to the impact of 
de-orbited past satellites (iCraddocld 120111 : iSchultz et al\ 
Il982t iBottke et al\ 12000(1 . Although such an estimate 
is limited because the elongated crater population could 
also include grazing asteroid impacts and may not ac¬ 
count for all prior martian satellites, the estimated mass 
of past and present satellites is still considerably greater 
than the current masses of Phobos and Deimos, which 
are 1.07 x 10^® and 1.48 x 10^® kg, respectively. 

Recently, the feasibility of forming the mar¬ 
tian moons from giant impacts has been studied 
both analytically (ICraddockI 12011(1 and numer ically 
([Rosenblatt and Charnozl 1201^ . ICraddocld (12011( 1 esti¬ 
mated that the mass of a circum-Mars debris disk should 
be twice the mass contained in the present and past mar¬ 
tian satellites {Mgats/Md ^ 0.5), based on the assump¬ 
tion that half of the disk material migrates inward while 
the other half migrates outward from Mars to form satel¬ 
lites. Comparing this estimate to the impactor mass of 
^ 0.02 MMars that would be necessary to cause the cur¬ 
rent martian spin-rate, Craddock estimated that ~ 0.01 - 
0.4 % of the impactor mass must be ejected into orbit to 
form the martian moons, a circum-Mars accretion disk of 
- Ixl0i®-5xl0i®kg {Md/Mimp -- IxlO-'i-dxlO"®). 

To determine if such an inserted mass could produce 
the total estimated mass of prior a nd cu rrent martian 
satellites, (Rosenblatt and Charnod (|2012ll modeled the 
evolution of a 10^® kg circum-Mars debris disk, a prob¬ 
lem significantly different from prior models o f satellite 
forma tion in a circum-Earth debris disk fe.g.. (Ida et al\ 
(1997(1 because of the l ow mass and density of t h e pro - 
posed martian disk. (Rosenblatt and Charnod (12012( 1 
modeled circum-Mars disk evolution in both a strong 
tidal regime, where accretion occurs at the Roche limit 
close to the planet, and a weak tidal regime, where 
accretion occurs furth er away from the plane t . In 
the strong tidal regime, (Rosenblatt and Charnozj (|2012[1 
found that to account for a moonlet population of 8.3 x 
10^^ — 1.5 X 10^® kg, the disk must have had a mass 
of 8.3 X 10^® — 1.5 X 1 0 ®^ kg , 30 — 80 times the mass 
predicted by iCraddocS (l201lh. If the relation, Md 
100 Msats, obtained bv iRoseiiblatt and Charno j (12012( 1 
is correct, in order to be consistent with Mi^p ^ 0.026 
MMars and Msats 1 X 10"® - 2 X 10“®, Md/M,mp 
must be (Msats/Mimp) x {Md/Msats) ^ 0.004 - 0.08. 
So far, the only other direct calculation of Md/Mimp is 
bv lCanup and Salmonl (|2014(1 . which was reported at the 
DPS meeting in 2014. 

To directly determine an estimate for Md/Mimp, we 
performed numerical Smoothed Particle Hydrodynamics 
(SPH) simulations of giant im pacts into Mars. We ta ke 
our canonical giant impact from lMarinova et al\ (12008( 1. a 
45 degree impact with an impact energy Eimp = 3 x 10®® 
J and Mimp = 0.026 MMars, which has been shown to 
produce a crustal dichotomy similar to that currently 
observed on Mars. To test the robustness of the result, 
we vary impactor angle and impact energy, particularly 
to account for other proposed mechanisms of dichotomy 


formation, such as an impact me gadome requiring an im¬ 
pact with Eimp ~ 1 — 3 X 10®® J (iReese et a/.ll20lill . Our 
numerical study of a range of potential impacts allows for 
improved inputs into circum-Mars disk evolution models. 

2. SPH SIMULATIONS 
2.1. SPH method 


We model giant impact scenario s onto Mars using 
the SPH method (jMonasfhanl (1992(1 . SPH codes have 
been used to model several types of giant i mpacts, in- 
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and other types of large-scale plairetary collisions 
(|Agnor and AsphanBl(2004t (Genda et a/.l(2012( l. An SPH 
code h as also been used to in odel the Borealis impact on 
Mars (|Marinova et al\ (2008(1 , but these simulations fo¬ 
cused on the structure of the impact basin and did not 
examine the formation of a debris disk. _ 

We use the SPH code described in (Genda et al\ ( 2012 1. 
which is based on the convention of (Monaghanl ( 1992 ). 
Our si mulations use t he Tillotson equation of state 
(EOS) (jTillotsonl (1962(1 . which is widely us ed in large- 
scale planetary iinpact simulations (e.g . . (Benz et al\ 
(1987(: (Canup et al\ (2001( : (Marinova et al\ (2008(1 . The 
Tillotson EOS parameters used in our study are identical 
to tho se used in prior simu lations of a giant impact on 
Mars (jMarinova et aL([2008(l . To represent Mars, we use 
a differentiated iron/mantle target body with a core to 
mantle mass ratio of 0.3 and a total weight of 6.0 xlO®® 
kg. The impactor is an undifferentiated basalt body with 
a mass of 0.026 x target mass. Our canonical impact has 
an impact velocity of 6 km/s, corresponding to an impact 
energ y of 3 x 10®® J, simila r to the nominal simulation 
from iMarinova et al\ (12008( 1. We vary impact velocity 
(a = Vimp/vesc) and impactor mass, Mimp, to obtain 
disk masses for several types of giant impacts into Mars. 

After each SPH simulation is complete, we quantify the 
mass and angular momentum of the impact-generated 
disk u sing the method described in iGnk and Steward 
(j2012ll : the disk contains SPH particles that have a neg¬ 
ative total energy and an equivalent circular orbital ra¬ 
dius that is greater than an estimated planetary ra¬ 
dius defined by a 1 g cm"® density contour. This 
method or a similar calculation is generally used to esti¬ 
mate disk masses from large collisions (e.g.,[Canu 3|2001 
iGnk and Stewartll2012l l. however, it is only an effective 
estimate for the initial mass of the disk shortly after im¬ 
pact. For timescales greater than several tens of hours af¬ 
ter impact, artificially induced spreading in the SPH sim¬ 
ulation leads to an ineffective disk mass calculation. Fur¬ 
thermore, the SPH simulation is limited in that the col- 
lisionless particle representation cannot accurately simu¬ 
late accretion in a disk, a computation better estimated 
using numerical models of disk evolu tion fe.g.. (Ida et al\ 
(1997(: [Rosenblatt and Charnozj(201^ . Therefore, our ini¬ 
tial disk mass estimates are meant to be compared with 
or used as inputs for long-term disk evolution models, 
and because the damping of inclination is rapid in com - 
parison to the disk’s radial evolution ((Ida et a/.l (1997(1 . 
disk mass should not change significantly beyond the 
time limit of our simulations due to damping. 

For our canonical impact, we conducted a resolution 
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TABLE 1 

Disk masses for several SPH simulations with various impact parameters^. 


Run 


N 

^imp (1^§) 

a 

0(°) 

E. 

imp (J) 

7 


Md 

(kg) ® 

Li 

0 

5 

X 

lO"' 

1.68 x 10^2 

1.4 

45 

3 

X 

1029 

0.028 

0.026 

4.3 

X 

1020 

0.95 

1 

1 

X 

105 

1.68 x 10^2 

1.4 

45 

3 

X 

1029 

0.028 

0.029 

4.8 

X 

1020 

0.98 

2 

3 

X 

10® 

1.68 x 10^2 

1.4 

45 

3 

X 

1029 

0.028 

0.033 

5.5 

X 

1020 

0.92 

3 

1 

X 

10® 

1.68 X 1022 

1.4 

45 

3 

X 

1029 

0.028 

0.033 

5.5 

X 

1020 

0.90 

4 

3 

X 

10® 

1.68 X 1022 

1.4 

30 

3 

X 

1029 

0.028 

0.024 

4.0 

X 

1020 

0.89 

5 

3 

X 

10® 

1.68 X 1022 

1.4 

60 

3 

X 

1029 

0.028 

0.017 

2.8 

X 

1020 

1.01 

6 

3 

X 

10® 

1.68 X 1022 

2 

45 

6.; 

! X 102® 

0.028 

0.024 

4.0 

X 

1020 

0.95 

7 

3 

X 

10® 

3.5 X 1022 

1.4 

45 

6 

X 

1029 

0.058 

0.033 

1.14 

X 

1021 

0.86 

8 

3 

X 

10® 

6.0 X 1022 

1.4 

45 

1 

X 

IqSO 

0.100 

0.023 

1.4 

X 

1021 

0.89 

9 

3 

X 

10® 

1.0 X 1023 

1.9 

45 

3 

X 

1030 

0.167 

0.013 

1.3 

X 

1021 

0.76 


^ Here, we list our ten simulations according to run number, number of particles in the simulation (N), 
impactor mass (Mimp), relative impact velocity (a = Vimp/vesc), impact angle (9), energy of the impact 
{Eimp), impactor-to-target mass ratio (7 = Mi^pIMtarget), disk-to-impactor mass ratio (Md/Mijjip), 
disk mass (M^), and normalized angular momentum (LJ). 

^ M^/Mtot, Md, and LJ are all reported at t = 30 hrs. 

^ We normalize the angular momentum as = Ld/(MdyJGMMarsRRoche)- 



hours 

Fig. 1.— Disk mass versus time for several total particle numbers. 

test (Runs 0 — 3, Table[T]) to compare disk mass evolution 
versus time for several total SPH particle counts (Fig. 
[T]). For the first 40 hours after the impact, the disk mass 
evolution of a simulation with 3 x 10® particles is similar 
to the simulation with 10® particles. Because of this, we 
use a resolution of N = 3 x 10® particles for the remainder 
of our simulations. 

2.2. SPH results 

Because the proposed Borealis impact is lower energy 
than the putative Moon-forming collision, impacts stud¬ 
ied in this investigation produced much smaller disks 
than previous studies of Earth-Moon formation. The 
disk masses for the collision scenarios we investigated 
are shown in Table [T] For our canonical impact (Run 2, 
Tabled]), a disk of mass 5.55 x 10^® kg is produced 30 
hours post-impact, composed of ^ 280 particles. 

An example of the evolution of a debris disk formed 
after our canonical Borealis-scale impact is shown in Fig. 
m In Fig. El the blue particles represent particles that are 
part of the planet or lack sufficient angular momentum 
to be considered part of the disk. The red particles are 
part of the disk, and black particles have been ejected 
from the system. The initial shock of the impact creates 
a large amount of ejecta that orbits the planet. While a 


degree of spreading of the ejecta occurs, a large amount 
of ejecta remains relatively close to the planet in a disk. 
An edge-on view of the disk is also shown in Fig. [3l The 
disk appears to undergo relatively quick damping, and 
Fig. dl shows that the disk mass becomes more constant 
as simulation time increases. 

To test the robustness of our result, we varied impact 
angle and impact energy, and compared the resulting de¬ 
bris disk masses. For impact angle, in addition to the 
nominal 45 degree impact angle, we also tested impacts 
with angles of 30 and 60 degrees. These changes in im¬ 
pact angle only slightly reduced the disk mass. For im¬ 
pact energy, we ran several simulations with higher im¬ 
pact energy by increasing the impactor mass and/or ve¬ 
locity (Runs 6—9, Tabled])- While increasing the impact 
energy did result in larger mass debris disks, the highest 
energy impact we examined (Run 9) only resulted in ap¬ 
proximately double the disk mass of the canonical impact 
(Run 2). In all cases, the amount of disk material pro¬ 
duced was about 1 to 4 % of the impactor mass. There¬ 
fore, we can conclude that Md/Mimp ^ (1 — 4) x 10“^ 
is robust for Borealis -scale impacts. Th i s app ears to 
agree with the result of lCanuo and SalmonI (|2014[ ). which 
found that an impactor with a few percent of Mars’ mass 
produced a disk several orders of magnitude more mas¬ 
sive than Phobos and Deimos. 

We also examined the disk angular momentum, scaled 
as L*J^ = Ld/ ( Mdy/GMMarsRRoche ) ■ If L*^ is greater than 
1, most of the mass is outside of the Roche limit, there¬ 
fore this quantity represents the compactness of the disk 
compared to the Roche limit. Our simulations (Table [IJ 
produce disks with < 1 for all but one case. We 
find that LJ decreases when both the impactor mass 
and impact energy increase. The higher angle 9 — 60° 
impact (Run 5) produces the least compact disk. Our 
canonical simulation (Run 2) produces L*^ = 0.92, repre¬ 
senting a slightly compact disk where most of the mass 
is within the Roche limit, but moderate mass is also 
present outside of the Roche limit. The presence of 
mass outside the Roche limit is also reported in the 
iCanuD and SalmonI (|2014[ 1 abstract. Because a value of 
TJ ^ 0.9 sh ould correspond to the s t rong tidal regime 
described bv iRosenblatt and Charnoi (I20I2I1 . our results 
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Fig. 2. — Snapshots of an SPH simulation during the impact 
and disk formation. Ejected particles are shown in black and disk 
particles are shown in red. The blue particles are considered part 
of the planet. 
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Fig. 3. — An edge on view of the simulation from Fig. 1 at f = 
4 and 44 hr. Colors of particles are the same as in Fig. 1. 

indicate that accretion in the strong tidal regime is fa¬ 
vored over accretion in the weak tidal regime. 

3. DISCUSSION AND CONCLUSION 

Our examination of a Borealis-scale impact (Run 2, 
Table[T]) shows that ^ 5 x 10^° kg of disk material is pro¬ 
duced 30 hours post-im pact, well within the 8.3 x 10^^ — 
1.5 X 10^^ kg disk mass IRosenblatt and CharnozI (I2012H 
require in their strong regime to produce sufficient satel¬ 
lite mass to explain the elongated crater distribution on 
Mars. Modest adjustments to the impact angle or veloc¬ 
ity produce similar disk masses, suggesting that sufficient 
material to produce Phobos and Deimos, or even the en¬ 
tire past/present satellite population, could have been 
ejected in a Borealis-scale impact. The larger impacts 
we examined, analogou s to collisions requ ired to form im¬ 
pact megadomes fe.g.. iReese et al.ll2011ll . produce over 
twice the disk mass as the smaller scale impacts. This is 
at the upper limit of the disk mass required to produce 
the theoretical moonlet population and greater than the 


disk mass necessary to produce Deimos in the weak ac¬ 
cretion regime. 

Our finding that 1 to 4 % of the impactor mass is in¬ 
serted in to an accretion disk is 20 — 100 times higher 
than the iCraddoctd (1201 Ih analytical estimate that 0.01 
to 0.2 % of the impactor i s inse rted into orbit. This 
is partly because iCraddocS (|2011ll estimates the ejected 
mass by assuming that half of a circum-Mars disk would 
migrate outward and accrete into satellites, and therefore 
the disk mass should be twice the mass of the past and 
present satellites. Additionally, the Gault and 7r-group 
scaling laws used to estimate the mass of past satellites 
are valid only for 90° impacts therefore should only pro¬ 
vide a lower limit on the masses of the sate llites that 
could have formed elongate craters on Mars (iCraddockI 
l2011h . By numeri cally modeling the evolution of c ircum- 
Mars debris disks. I^senblatt and Charnod (|2012fl found 
that only 1% of the disk mass necessarily accreted into 
satellites. Therefore, although o ur simulati o ns pr oduce 
disks with higher mass than the iCraddocS (l2011f l esti¬ 
mates, if less of the disk mass accretes into moonlets and 
satellites, then similar satellite mass distributions can be 
achieved. 

It should be noted that our study examines a narrow 
parameter space of Borealis-scale impacts that we find 
produce di sks more massive than the 10^ ^ kg disk ex¬ 
amined by IRosenblatt and Charnod (l2012fl . Therefore, 
altho ugh we assume that 1% of the di sk mass forms satel¬ 
lites (|R,osenblatt and CharnozI 12012( 1. it is not certain 
that this estimate holds at larger disk masses, and further 
simulations of circum-Mars disk evolution for larger disk 
masses are necessary. Additionally, accretion in more 
massive disks (> 10^® kg) coul d produce satellites larger 
in siz e than Phobos or Deimos (|R,osenblatt and CharnozI 
unii); however, these satellites could eventually de-orbit 
and impact Mars, depending on where they form in the 
disk. Furthermo re, in addition to Borealis-scale impacts, 
iCraddocS (|201ir ) suggests that the Elysium and Utopia 
basins could have been formed by large enough impacts 
to produce the current martian spin rate, suggesting im¬ 
pacts as low as 0.015 Mars masses. A basic power law ht 
to our simulations with fixed a = 1.4 and 0 = 45 (Runs 
2, 7, and 8) yields Mj, = 8 x 10^^(Mimp)°'^^, meaning 
an impactor with 0.015 target masses would produce a 
disk of mass ~ 3 x 10^° kg. Direct simulations of smaller 
impacts would be useful, but as the impact energy de¬ 
creases, more total particles are needed in the simulation 
if the final disk is to contain at least the couple hundred 
particles necessary to accurately estimate the disk mass. 
Because our impact en e rgies fall in the most-likely range 
reported bv ICraddockI (|2(111[ 1. we leave the examination 
of less energetic impacts for future work. 

Our simulations show that for Borealis-scale impacts, 
enough material is ejected into orbit to form accre¬ 
tion disks that could produce martian satellites. Our 
result is inconsistent with the formation of only Pho¬ 
bos and Deimos in the weak tid al regime considered 
by IRosenblatt and Charnod (|2012f) because our disk is 
100 — 1000 times more massive and the values of our 
disks mostly indicate too much compactness for the weak 
tidal regime. Therefore, our results support accretion in 
the strong tidal regime. Although additional simulations 
of disk evolution with disks of initially larger mass are re- 
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quired to further investigate scenarios of martian moon 
formation, our results provide a numerical estimate of 
disk mass production in Borealis-scale impacts. 
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